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We develop a method based on the relativistic coupled-cluster theory to incorporate a perturbative 
interaction to the no-pair Dirac-Coulomb atomic Hamiltonian. The method is general and suitable 
to incorporate any perturbation Hamiltonian in a many electron atom or ion. Using this perturbed 
relativistic coupled-cluster (PRCC) theory we calculate the electric dipole polarizability, a, of Neon. 
The linearized PRCC results are in very good agreement with the experimental value. However, the 
results of the nonlinear PRCC shows larger uncertainty but it is consistent with the observations 
from earlier works. 
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I. INTRODUCTION 

The electric dipole polarizability, a, of an atom or 
molecule is an important property as it describes the in- 
teraction with an external electric field. Knowing accu- 
rate value of atomic a is imperative in different areas of 
research. To mention a few, the parity nonconservation 
in atoms [l| and the ultracold atoms are of cur- 

rent interest. In short, knowledge of a is essential in any 
experiment which involves the atom interacting with an 
external electric field. Theoretically, a is calculated using 
the first order time- independent perturbation theory. A 
recent review Q gives a detailed account of the current 
status of the theoretical results of atomic and ionic po- 
larizabilities. An excellent source of information related 
to atomic polarizability is ref. [6[|. 

In the present work we extend the standard relativis- 
tic coupled-cluster (RCC) theory to include an additional 
perturbation Hamiltonian. Here, it must emphasized 
that the coupled-cluster (CC) theory 0, H| is one of the 
most reliable quantum many body theory. We refer to 
the CC theory with the additional perturbation as PRCC 
method. The theory is flexible and formulated to incor- 
porate multiple perturbations which may have different 
tensor structures in the electronic sector. The perturba- 
tion could be internal, like hypcrfine interaction or exter- 
nal like the static electric field. For the present work, as 
a is the quantity of interest, we take an external electric 
field as the perturbation. The coupled-cluster theory has 
been widely used for atomic I9l. [l0ll. molecular 11 , nuclear 
[lU, condensed matter physics |13f calculations. The pre- 
vious approaches of CC theory for atomic calculations are 
conceptually different from the PRCC method in one as- 
pect. In PRCC, the CC single and double (CCSD) excita- 
tion operators within the electron sector could be tensor 
operators of higher ranks. In which case the associated 
angular factors are complicated and are evaluated using 
diagrammatic method [lj]. The PRCC equations are, 
at first glance, complicated in appearance as it involves 
two sets of cluster operators. However, the equations are 
linear and not nonlinear as in the full RCC theory. For 



calculations, we used the equivalent of standard CC sin- 
gle and double but the theory can be extended to higher 
orders. 

The paper is organized as follows. In the Section. II 
we provide a brief description of the PRCC method and 
examine the structure of the cluster equations. This is 
followed with a short presentation of the details of the 
numerical aspects and actual models used in Section. III. 
The Section. IV discusses the results and we end with 
the conclusions. 



II. THEORETICAL METHODS 

The Dirac-Coulomb Hamiltonian H uc , consisting of 
relativistic one-electron terms and electrostatic electron- 
electron interaction, is the starting point of our calcula- 
tion. For an iV-electron atom 

N 1 

i? DC = 5>a, • Pj + (A - i)c 2 - v N ( n )] + (i) 

1—1 Z<J J 

where, a^s and the fii are the Dirac matrices, V/v(?"i) is 
the nuclear potential and the last term represents the 
electron-electron Coulomb interactions. For a closed- 
shell atom the eigen- value equation is 



(2) 



In the RCC theory the ground state wave-function of an 
atom is written as 



l*o)= e T( >o>, 



(3) 



here |$ ) is the ground state reference state. In the 
CCSD approximation the cluster operator T' ' = T^ ' + 

T 2 . The advantage of CC theory is, it is an all or- 
der non-perturbative theory and ensures that electron- 
electron correlation diagrams is accounted fully. For fu- 
ture reference T^°> are referred to as the unperturbed 
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cluster operator. In the CCSD approximation the ampli- 
tude equations are written as 

(4) 
(5) 



&S|ff N |$o> = 0, 
(*S|F N |$o> = 0. 



Here = e~ T< ' H-^e Ti0) is the similarity transformed 
Hamiltonian and H N = H BC - (<& \H DG \$ ) is the nor- 
mal ordered Hamiltonian fl5j . 

In the presence of an additional interaction, which 
could be an internal or external perturbation Hamilto- 
nian H' , \^&o) is modified. Examples of internal per- 
turbations are hyper-fine interaction, Breit interaction, 
etc. And examples of external perturbations are Stark 
effect, Zeeman effect, etc. The perturbation Hamilto- 
nian, in the electron sector, can be a tensor operator of 
rank one or higher. The effect of these perturbations, al- 
though important, are less in magnitude than the resid- 
ual Coulomb interaction and a first order treatment is 
sufficient. To incorporate H', we have developed a per- 
turbed RCC (PRCC) theory. 

For the present work we take an external electric field 
as the perturbation in an N-clectron atom. The atomic 
Hamiltonian is then modified to 



where, AB represents all possible contractions between 
the two operators A and B. These are non- linear cou- 
pled equations. However, are solutions of the unper- 
turbed RCC equations [17|, which are solved separately 
With this consideration, the T' 1 ) equations are reduce to 
a set of linear algebraic equations. The terms on the left 
hand side of Eq. ([5]) and © can be simplified to 



(0)rA(l)l 



{jffNT (0) T (0) T (i) } = {HnT W t W t w } 



where, {A} represents the normal order form of the op- 
erator A. Similarly, the terms on the right hand side 
expand to 



(0)^(0)^,(1) 



{H int T^} = {ffintTD + {HintT?'}, 



'(0)i 



H = H 



DC 



(6) 



where Hi nt = ^. i=1 .E, with E as the external elec- 
tric field and A is the perturbation parameter. Eigen- 
value equation with the modified Hamiltonian is H \^o) = 
Ei\^o). In the PRCC theory the new ground state wave- 
function is 



{HintT (o) T (o) } = {HintT p T (0) } + {H . ntT i0) T (0 )} _ 

For further evaluation we use diagrammatic techniques 
to calculate the contributions from these terms. In total 
there are 42 single diagrams and f 02 doubles diagrams. 

From the time independent perturbation theory the 
dipole polarizability of an atom is 



a = -2 



l*o> 



T (0) +AT (1) 



*0>- 



I 



OP |D|tt 7 )(tt f |D|W ) 
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Here we have introduced the PRCC operator T 1 - 1 ), which 
incorporates the many-body effects of H int . It operates 
on the electronic sector and is a rank one operator as the 
electronic part of H lnt is a tensor operator of rank one. 

To first order in A we have |^o) = 
in the CCSD approximation : 



where, \^fj) are the intermediate atomic states. In the 
PRCC theory the sum over states is implicit and a 

is the expectation value of the operator e T< ' De T< ) + 

e T<0) _De T< ' . After expansion, collecting terms upto scc- 



l_l_ i yp(i)j|\|/ \ anc [ ono - or dcr in cluster operators 



,(i) 



T^ ; . More 



detailed description on the structures of and T^ 1 ' 
operators are given our previous publication [16j . 
The PRCC amplitude are solutions of the equations 



= ^($o|{TPb} + {Dlf >} + {T^BT^} 



(*SI 
1 
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.{T a C°)t DT W} + {T^-DTP} + {Ti^tDT^} 
-{T 2 MiYti 0) } + {Ti^tDTW} + { T 2 (DtD T (0) } 

.{r 2 (o)tDf«}|$ }, (ii) 



fl N T'(o)T'(o)j'(o) T (i) 



l*o> = -(*SI 



2! 



H int T^T 



(0)t(0) 



l*0>, 



(8) 



21' 



l*o>=-<*SI 



ffmtT<°> 



l*o>, 



(9) 



where N = ($ |{Ti (0)t Ti(°)} + {T 2 ^T 2 ^}\^ ) is the 
normalization constant. The unperturbed wave-function 
is normalized, but we have to normalize the perturbed 
wave- function. 



III. CALCULATIONAL DETAILS 

Using PRCC we have done a systematic study of dipole 
polarizability of Neon. For precision atomic theory cal- 
culation the orbital basis set is an important factor. In 
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our calculations we use even tempered Gaussian type or- 
bitals (GTOs) [l8j]. The radial part of the Dirac bi-spinor 
is expressed in terms of the linear combination of GTOs. 
The GTOs 



r'" n e 



(12) 



The index p is the number of basis functions. The ex- 
ponent a p depends on two parameters ao and j3 and are 
related as a p = ao/3 p_1 , where p = 0, 1 . . . m and m is the 
number of g^ p (r) considered. For present calculations we 
use even tempered basis, where ao and /3 are unique for 
each of the symmetries. The symmetry wise values of the 
parameters are listed in Table. |TJ Although, in principle, 



TABLE I. The GTO even tempered basis parameters used in 
the present calculations. 



Symmetry 

P 

Symmetry 
ao 



Sl/2 Pl/2 P3/2 ^3/2 d 5 / 2 fz,/2 

0.0925 0.1951 0.1917 0.0070 0.0070 0.0069 

1.45 2.71 2.71 2.70 2.70 2.69 

J7/2 #7/2 59/2 

0.0069 0.0069 0.0069 
2.69 2.69 2.69 



a complete set of orbitals are needed, it is near impossi- 
ble to go beyond a few hundred. Even at a few hundred 
the computational requirements is very high. Another 
practical consideration, with further increase the gain in 
accuracy is marginal or non-existent once the basis set 
converges. The basis parameters are optimized such that 
the core orbital energies are in good agreement with the 
GRASP92 [l9| results. For information, the orbital en- 
ergies are listed in Table. HU 



TABLE II. Orbital energies of obtained from GRASP92 and 
GTO in atomic units. 



Orbital 


lSl/2 


2S1/2 


2pi/ 2 


2f>3/2 


GTO 


-32.8177 


-1.9357 


-0.8526 


-0.8480 


GRASP92 


-32.8145 


-1.9387 


-0.8528 


-0.8482 



To test and check the theory, we consider the linearized 
PRCC theory. There are then 10 singles and 10 doubles 
diagrams each. However, only 6 of the singles diagram 
but all the doubles diagrams contribute when DHF or- 
bitals are used. Detailed descriptions of the diagram- 
matic calculations are given in rcf. (l6j . The result, 
along with previous and experimental values, are given 
in Table. IIII1 It shows that our results agrees very well 
with the experimental data and indicates that the PRCC 
theory, even at the linear level, gives accurate results for 
a single reference system like Ne. To analyze the impact 



TABLE III. Contribution from Linearized PRCC to the static 
dipole polarizability of Ne and comparison with previous re- 
sults. 

This work CCSDT[2CV] RCCSDT[2TJ MBPT[22J Expt.[23j 
2.6695 2.6648 2.697 2.665 2.670(5) 



of basis set truncation, we examine the convergence of 
a with the size of basis set. For this we start with a 
basis set of 50 GTOs and do a scries of calculations by 
increasing the basis size in steps. The value of a con- 
verges to 2.6695 when the basis set size is 124. However, 
for confirmation we increase the basis set size upto 171 
and results are listed in Table. IIVI In this calculation we 



TABLE IV. Convergence pattern of a of Ne as a function of 
the Basis set size. 



No. of orbitals 


Basis size 


Polarizability 


50 


(10s,6p,6d,4/, 4fl) 


2.7279 


60 


(12s,7p,7d,5/,5 9 ) 


2.7087 


75 


(13s, 9p, 9d, 7f, &g) 


2.6849 


91 


(15s, lip, lid, 8/, 85) 


2.6712 


108 


(20s, 13p, lid, 11/, 9g) 


2.6696 


124 


(22s, Up, Ud, 13/, 10fl) 


2.6695 


145 


(27s, 17p, 16d, 14/, 12 5 ) 


2.6695 


163 


(29s, 21p, Yld, 16 f, 13g) 


2.6695 


171 


(31s, 23p, 18d, 16/, 13g) 


2.6695 



have considered finite size Fermi density distribution for 
the nucleus. 



=(r) 



Pa 



1 + e O-c)/a ' 



(13) 



Here, a = i41n(3). The parameter c is the half charge 
radius so that p n uc(c) = po/2 and t is the skin thickness. 
The PRCC equations are solved iteratively using Jacobi 
method, we have chosen this method as it is paralleliz- 
able. The method, however, is slow to converge. So, we 
use direct inversion in the iterated subspace (DIIS)[24[ 
to accelerate the convergence. 



IV. RESULTS AND DISCUSSIONS 



In the properties calculations the CC expression of 



d)t 



the polarizability operator, e De + 



n(0)t 



De 1 



is a nontcrminating series. However, as described ear- 
lier, in the present calculations we consider upto second 
order in T. The contributions from the higher order 
terms, based on previous studies with an iterative all 
order method[25l|. is negligible. The contributions from 
different terms in Eq. (fTTj) are listed in Table. [V] As 
evident from the table, the dominant contribution arises 

from {T\ D} and its hermitian conjugate. This is not 
surprising as these terms subsume the DF contribution 
and core-polarization effects. The general trend is, for 
closed-shell atoms, the DF and core-polarization effects 
are the leading order and next to leading order, respec- 
tively. Coming to the pair correlation effects, the leading 

contribution arise from {T 2 (1)t DT 2 (0) } and its hermitian 

conjugate. This is along the expected lines as the 

amplitude is larger, compared to T^°\ on account of pair- 
correlations. The contributions from the remaining terms 
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TABLE V. Contribution to a of Ne from different terms of 
the dressed dipole operator in the linearized PRCC theory 



Contributions from a 

{ T^ 1)t b} + h.c. 2.6610 

{ Ti (1)t br 2 (Q) } + h.c. -0.0478 

{ Ti (1)t bf 1 ( ° 1 } + h.c. 0.0644 

{ T 2 (1)t E)f 1 ( " ) } + h.c. -0.0062 

{T 2 (1)t bf 2 (0) } + h.c. 0.0961 

Normalization 1.0367 

Total 2.6695 



are small and cancellations reduce the combined contri- 
bution even further. 

The next level of calculation is to consider all the terms 
in the Eq. and ([§]) , which we refer to as the non-linear 
PRCC. The term wise contributions are listed in Table. 

TABLE VI. Contribution from all the terms in the nonlinear 
PRCC theory. 



CC terms a 

{ T[ in i>} + h.c. 2.7344 

{ T 1 (1)t br 2 (Q) } + h.c. -0.0492 

{Ti (1)t E)f 1 (0) } + h.c. 0.0670 

{ T a (1)t br i (0) } + h.c. -0.0058 

{T 2 (1)t bf 2 (0) } + h.c. 0.0924 

Normalization 1.0367 

Total 2.7383 



IVII and the net result of 2.7383 is 2.6% larger than the 
linearized PRCC result. As evident from the table, most 
of the change is attributed to {T\ D} and hermitian 
conjugate. Contribution from this term is 2.7% larger in 
the nonlinear PRCC, which is comparable to the change 
in the value of a. This is one of the case where higher 
order calculations does not translate into improved ac- 
curacy. A similar situation, but in a different context, 
was observed in a detailed analysis of contributions from 
nonlinear terms in the CCSD and dressing to calculate 
the magnetic dipole hyperfine constant of Li J26|. Like 
in the work referred, the contributions from higher order 
cluster operators, triple and quadruple excitations, could 
be of different phase and bring a closer to experimental 
data. However, to examine or confirm this requires more 
detailed analysis and is beyond the scope of the present 
work. 

Through a series of rigorous calculations, we examine 
the changes in a, and associate it with a nonlinear term in 
Eq. ([SJ and ©. At the second order, there is an anoma- 
lously large contribution from {H-^T^T^}, it induces 
a changes of 0.0808 a.u. to the net result of a. This term 



TABLE VII. Two of the leading order terms in the nonlinear 
PRCC Theory. 



Contributions {H N T 2 (0) T^} {HnT^T^T^} 



{T< in b} + h.c 


2.7456 


2.6628 


{TiWtDT 2 (0) } + h.c. 


-0.0492 


-0.0478 


{T 1 ( 1 )tbT 1 (0) } + h.c. 


0.0674 


0.0642 


{T 2 WtbT 1 (0) } + h.c. 


-0.0058 


-0.0058 


{T 2 (1)t br 2 (0) } + h.c. 

Normalization 
Total 


0.0933 
1.0367 
2.7503 


0.0922 
1.0367 
2.6677 



accounts for the large change of a in the nonlinear PRCC 
calculations. Compared to this term, the contribution 
from other terms at this order are marginal. The next 

largest contribution arises from {HnT^T^}, it con- 
tributes 0.0086 a.u. Th e other contributio ns arc 0. 0004 

and 0.0034 a.u. from {H^jf^T^} and {H N T^ 0) T^}, 
respectively. 

At the third order {H^T^T^T^} and 

{H^T^T^T^} contribute equally, 0.0077 a.u. 
each. The contribution from the last term at this order, 

{HnT^T^T^}, is -0.0018 a.u. To illustrate the 
relative changes arising from the third order terms, we 
list the contributions from the leading order terms in the 
second order and third order in Table. IVIII It is evident 
from the table that the difference between the second 
and third order contributions arises from the {T^D} 
and its hermitian conjugate. 

At the fourth order there is only one term 

{HnT^T^T^T^} and contributes 0.0077 a.u. This 
detailed study implies that the higher order terms in the 
PRCC equations, third and fourth order, have negligible 
effect on the electric dipole polarizability. Since the effect 
of the higher terms are tightly coupled to the electron 
correlation effects, a similar trend may occur in other 
properties as well. 

To estimate the uncertainty in our calculations, we 
have identified two sources in the calculations using 
PRCC with CCSD approximation. First type of error 
is associated with the orbital basis set truncation and 
the termination of iteration while solving the cluster am- 
plitudes. Based on the basis set convergence, as de- 
scribed earlier, the uncertainty from the basis set trun- 
cation is negligible. Similarly, the uncertainty from the 
termination of cluster amplitude calculation is negligi- 
ble as we set 10 -6 as the convergence criterion. The 
second type of error arises from the truncation of the 
CC theory at double excitations and the truncation of 
e T<1) De Tm + e T<0> De T<1) . Based on other detailed 
studies on the contributions from the triples and quadru- 
ple excitations could be in the range of «-2.6%. So that 
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TABLE VIII. The contribution to a from the fourth order 
term in nonlinear PRCC theory. 



Contributions From {HnT^T^T^T^} 



{T^Dj+h.c. 


2.6688 


{T^Wbff'j + h.c. 


-0.0478 


{T 1 ( 1 )tDT 1 W} + h.c 


0.0645 


{TaWtbT^^ + h-c 


-0.0062 


{T 2 (1)t br 2 (0) } + h.c 

Normalization 
Total 


0.0962 
1.0367 
2.6772 



it balances the larger error arising from the inclusion of 
the nonlinear terms in the PRCC theory. Again, based 
on other studies with iterative method 25| to incorporate 
higher order terms in the properties calculations with CC 
theory, the contributions from the third or higher order 



T' 1 ' 1 7-. T<°> 

me L>e 



e T<0) De T(1) is negligibly small. The 



contribution from Breit and QED corrections could be 
another source of error. However, as Za <C 1, where a is 
the fine structure constant for Ne the uncertainty from 
excluding Breit and QED correction could easily be sub 
0.01%. This is consistent with the estimates of the con- 
tribution from the Breit interaction to correlation energy 
[27I ] . Combining all the sources of error, the uncertainty 
for the calculations with nonlinear PRCC is w 2.6%. But 
due to possible fortuitous cancellations, the uncertainty 



with linearized PRCC calculations is sub 0.1%. 



V. CONCLUSION 

The PRCC theory provides a consistent approach 
within CC theory to calculate atomic properties. Al- 
though, in this paper we have demonstrated the use of 
PRCC to calculate a of Nc. The method is applicable 
to any atomic property. In PRCC theory, the number 
of cluster amplitudes is larger than the RCC theory but 
the PRCC equations are linear. The main feature of the 
method is accounting of all possible intermediate states 
within the basis set chosen for the calculations. With 
converged basis sets, the intermediate states incorporated 
in the properties calculations is practically complete. 

From the detailed calculations presented in this work, 
it can be concluded that inclusion of the nonlinear terms 
does not improve the accuracy of the properties calcu- 
lated with PRCC theory. However, The linearized PRCC 
provides results which are in excellent agreement with the 
experimental data. The lower accuracy with the nonlin- 
ear PRCC may be attributed to large cancellations be- 
tween contributions from different nonlinear term s. As 
to be expected, the leading order term is {Tp^D}. It 
accounts for, in the linearized PRCC calculations, 96% 
of the total value. 
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